A new non-linear least squares algorithm for the seismic inversion problem
نویسندگان
چکیده
We describe the problem of simultaneous determination of hypocentre locations and velocity structure, using data from microearthquake networks, and briefly discuss some of the techniques which have been proposed for its solution. In particular we examine the method of Thurber, as implemented in his program SIMUWL, where the problem is formulated as a non-linear least-squares problem. A new program SINEW is proposed which is a modification of SIMUL3L in which a somewhat different algorithm is used t o solve the underlying least squares problem. At each iteration, the vector of parameter adjustments is calculated inexactly via an iterative scheme which avoids explicit factorization of the (sparse) Jacobian matrix.
منابع مشابه
Exact and approximate solutions of fuzzy LR linear systems: New algorithms using a least squares model and the ABS approach
We present a methodology for characterization and an approach for computing the solutions of fuzzy linear systems with LR fuzzy variables. As solutions, notions of exact and approximate solutions are considered. We transform the fuzzy linear system into a corresponding linear crisp system and a constrained least squares problem. If the corresponding crisp system is incompatible, then the fuzzy ...
متن کاملNon-linear least-squares inversion with data-driven Bayesian regularization
The non-linear inverse problem is formulated in a Bayesian framework. The multivariate normal distribution is assumed in both the noise and prior distributions. However, only the structures of the covariance matrices have to be specified, estimation of the variance levels is included in the inversion procedure. The maximum a posteriori approximation is derived, and the final result is a weighte...
متن کاملNon-linear stochastic inversion of 2D gravity data using evolution strategy (ES)
In the current work, a 2D non-linear inverse problem of gravity data is solved using the evolution strategies (ES) to find the thickness of a sedimentary layer in a deep-water situation where a thick sedimentary layer usually exists. Such problems are widely encountered in the early stages of petroleum explorations where potential field data are used to find an initial estimate of the basin geo...
متن کاملLarge-scale Inversion of Magnetic Data Using Golub-Kahan Bidiagonalization with Truncated Generalized Cross Validation for Regularization Parameter Estimation
In this paper a fast method for large-scale sparse inversion of magnetic data is considered. The L1-norm stabilizer is used to generate models with sharp and distinct interfaces. To deal with the non-linearity introduced by the L1-norm, a model-space iteratively reweighted least squares algorithm is used. The original model matrix is factorized using the Golub-Kahan bidiagonalization that proje...
متن کاملA new stochastic 3D seismic inversion using direct sequential simulation and co-simulation in a genetic algorithm framework
Stochastic seismic inversion is a family of inversion algorithms in which the inverse solution was carried out using geostatistical simulation. In this work, a new 3D stochastic seismic inversion was developed in the MATLAB programming software. The proposed inversion algorithm is an iterative procedure that uses the principle of cross-over genetic algorithms as the global optimization techniqu...
متن کامل